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The Landau-Lifshitz-Gilbert (LLG) equation is a fascinating nonhnear evolu- 
tion equation both from mathematical and physical points of view. It is related to 
the dynamics of several important physical systems such as ferromagnets, vortex 
I filaments, moving space curves, etc. and has intimate connections with many of the 

well known integrable soliton equations, including nonlinear Schrodinger and sine- 
^ I Gordon equations. It can admit very many dynamical structures including spin 

waves, elliptic function waves, solitons, dromions, vortices, spatio-temporal pat- 
terns, chaos, etc. depending on the physical and spin dimensions and the nature of 
interactions. An exciting recent development is that the spin torque effect in nano- 
ferromagnets is described by a generalization of the LLG equation which forms a 
basic dynamical equation in the field of spintronics. This article will briefly review 
these developments as a tribute to Robin BuUough who was a great admirer of the 
' LLG equation. 
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Spin systems generally refer to ordered magnetic systems. Specifically, spin angu- 



^ , 1. Introduction 

in 

I lar momentum or spin is an intrinsic property associated with quantum particles, 

which does not have a classical counterpart. Macroscopically, all substances are 
magnetic to some extent and every material when placed in a magnetic field ac- 
quires a magnetic moment or magnetization. In analogy with the relation between 
the dipole moment of a current loop in a magnetic field and orbital angular momen- 
tum of a moving electron, one can relate the magnetic moment/magnetization with 
the expectation value of the spin angular momentum operator, which one may call 
simply as spin. In ferromagnetic materials, the moment of each atom and even the 
average is not zero. These materials are normally made up of domains, which ex- 
■ hibit long range ordering that causes the spins of the atomic ions to line up parallel 

to each other in a domain. The underlying interaction (Hillebrands & Ounadjela 
2002) originates from a spin-spin exchange interaction that is caused by the overlap- 
ping of electronic wave functions. Additional interactions which can influence the 
magnetic structures include magnetocrystalline anisotropy, applied magnetic field, 
demagnetization field, biquadratic exchange and other weak interactions. Based on 
phenomenological grounds, by including effectively the above type of interactions. 
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Landau and Lifshitz (1935) introduced the basic dynamical equation for magneti- 
zation or spin S{f, t) in bulk materials, where the effect of relativistic interactions 
were also included as a damping term. In 1954, Gilbert (2004) introduced a more 
convincing form for the damping term, based on a Lagrangian approach, and the 
combined form is now called the Landau-Lifshitz-Gilbert (LLG) equation, which is 
a fundamental dynamical system in applied magnetism (Hillebrands & Ounadjela 
2002; Mattis 1988; Stiles & Miltat 2006). 

The LLG equation for the unit spin vector S{f,t), because of the constancy 
of length, is a highly nonlinear partial differential equation in its original form for 
bulk materials. Depending on the nature of the spatial dimensions and interactions, 
it can exhibit a very large variety of nonlinear structures such as spin waves, el- 
liptic function waves, solitary waves, solitons, lumps, dromions, bifurcations and 
chaos, spatiotermporal patterns, etc. It exhibits very interesting differential geo- 
metric properties and has close connections with many integrable soliton and other 
systems, for special types of interactions. In the general situations, the system is 
highly complex and nonintcgrablc. Both from physical and mathematical points of 
view its analysis is highly challenging but rewarding. 

One can also deduce the LLG equation starting from a lattice spin Hamiltonian, 
by postulating appropriate Poisson brackets, and writing down the corresponding 
evolution equations and then introducing the Gilbert damping term phenomenolog- 
ically. Thus we can have LLG equation for a single spin, a lattice of spins and then 
the continuum limit in the form of nonlinear ordinary differential equation(ODE), 
a system of coupled nonlinear ODEs and a nonlinear partial differential equation, 
respectively, for the unit spin vector(s). Analysis of the LLG equation for discrete 
spin systems turns out to be even harder than the continuum limit due to the na- 
ture of nonlinearity. However, apart from exact analytic structures, one can also 
realize the onset of bifurcations, chaos and patterns more easily in discrete cases. 
Thus the LLG equation turns out to be an all encompassing nonlinear dynamical 
system. 

The LLG equation has also close relationship with several other physical sys- 
tems, for example motion of a vortex filament, motion of curves and surfaces, cr- 
models in particle physics, etc. One of the most exciting recent developments is that 
a simple generalization of the LLG equation also forms the basis of the so called 
spin torque effect in nanoferromagnets in the field of spintronics. 

With the above developments in mind, in this article we try to present a brief 
overviciw of the; different aspects of the LLG equation, concentrating on its nonlinear 
dynamics. Obviously the range of LLG equation is too large and it is too difficult 
to cover all aspects of it in a brief article and so the presentation will be more 
subjective. The structure of the paper will be as follows. In Sec. 2, starting from the 
dynamics of the single spin, extension is made to a lattice of spins and continuum 
systems to obtain the LLG equation in all the cases. In Sec. 3, we briefly point out the 
spin torque effect. Sec. 4, deals with exact solutions of discrete spin systems, while 
Sec. 5 deals with continuum spin systems in (1-t-l) dimensions and magnetic soliton 
solutions. Sec. 5 deals with (2-1-1) dimensional continuum spin systems. Concluding 
remarks are made in Sec. 6. 
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2. Macroscopic Dynamics of Spin Systems and LLG 

Equation 

To start with, in this section we will present a brief account of the phenomenological 
derivation of the Landau-Lifshitz-Gilbert (LLG) equation starting with the equa- 
tion of motion of a magnetization vector in the presence of an applied magnetic field 
(Hillebrands & Ounadjela, 2002). Then this analysis is extended to the case of a lat- 
tice of spins and its continuum limit, including the addition of a phenomenological 
damping term, to obtain the LLG equation. 



(a) Single Spin Dynamics 

Consider the dynamics of the spin angular momentum operator 5 of a free 
electron under the action of a time-dependent external magnetic field with the 
Zeeman term given by the Hamiltonian 

Hs = S.B{t), B{t) = f,oH{t), (1) 

where g, fiB and /xq are the gyromagnetic ratio, Bohr magneton and permeability in 
vacuum, respectively. Then, from the Schrodinger equation, the expectation value 
of the spin operator can be easily shown to satisfy the dynamical equation, using 
the angular momentum commutation relations, as 

1 < s{t) > = ^ < S{t) > xBit). (2) 

Now let us consider the relation between the classical angular momentum L of a 
moving electron and the dipole moment M.^ of a current loop immersed in a uniform 
magnetic field, M.^ = 2m ^' where e is the charge and m is the mass of the electron. 
Analogously one can define the magnetization M. = < 5>=7 < S >, where 
7 = ^j^- Then considering the magnetization per unit volume, M, from (2) one 
can write the evolution equation for the magnetization as 

^ = -^o[M{t)xH{t)], (3) 

where B = ji^H and 70 = /^oT- 

Prom (3) it is obvious that M.M = constant and M.H = constant. Conse- 
quently the magnitude of the magnetization vector remains constant in time, while 
it precesses around the magnetic field H making a constant angle with it. Defining 
the unit magnetization vector 

S(t) = S^ = l, S = (S'=,Sy,S'), (4) 

\M{t)\ 

which we will call simply as spin hereafter, one can write down the spin equation 
of motion as (Hillebrands & Ounadjela 2002) 

^ = -^o[S{t)xH{t)], H = {H\Hy,H^) (5) 
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Figure 1. Evolution of a single spin (a) in the presence of a magnetic field and (b) when 

damping is included. 

and the evolution of the spin can be schematically represented as in Fig. 1(a). 

It is well known that experimental hysteresis curves of ferromagnetic substances 
clearly show that beyond certain critical values of the applied magnetic field, the 
magnetization saturates, becomes uniform and aligns parallel to the magnetic field. 
In order to incorporate this experimental fact, from phenomenological grounds one 
can add a damping term suggested by Gilbert (2004) so that the equation of motion 
can be written as 



^= -7o[5(i)x^(t)]+A7o 



S X — 
dt 



, A <C l{damping parameter). (6) 



On substituting the expression (6) again for ^ in the third term of (6), it can be 
rewritten as 

(1 + A^o)^ = -jo[SxH{t)]-X-foSx[SxH{t)]. 
After a suitable rescaling oft, Eq.(6) can be rewritten as 

^ = (5 X H{t)) + XSx[Sx H{t)] 

= SxH^ff. (7) 
Here the effective field including damping is 

Heff = Hit) + X[SxH]. (8) 

The effect of damping is shown in Fig 1(b). Note that in Eq.(7) again the constancy 

of the length of the spin is maintained. Eq.(7) is the simplest form of Landau- 
Lifshitz-Gilbert equation, which represents the dynamics of a single spin in the 
presence of an applied magnetic field H{t). 

(6) Dynamics of LaMice of Spins and Continuum, Case 

The above phenomenological analysis can be easily extended to a lattice of 
spins representing a ferromagnetic material. For simplicity, considering an one di- 
mensional lattice of N spins with nearest neighbour interactions, onsite anisotropy. 
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demagnetizing field, applied magnetic field, etc., the dynamics of the ith spin can 
be written down in analogy with the single spin as the LLG equation, 

^ = SiXH.ff, i = l,2, N (9) 

where 

Heff = {S,+i + + AS^n^ + BSfHy + CS^n, + Hit) + ) 

- A{5,+i + Si-i + ASffi^ + BSffiy + CSfn, + H{t) + } x Si (10) 

Here A,B,C are anisotropy parameters, and nx,ny^nz are unit vectors along the 
X, y and z directions, respectively. One can include other types of interactions like 
biquadratic exchange, spin phonon coupling, dipolc interactions, etc. Also Eq.(9) 
can be generalized to the case of square and cubic lattices as well, where the index 
i has to be replaced by the appropriate lattice vector i. 

In the long wavelength and low temperature limit, that is in the continuum 
limit, one can write 

S^{t) = S{r,t), f= {x,y,z), 



a 



2 



Si+i + = S{r, t) + a.yS + y V^S' + higher orders, (11) 

(o: lattice vector) so that the LLG equation takes the form of a vector nonlinear 
partial differential equation (as a — >^ 0), 



dS{r,t) 

at 



--S X [{v^S + A5^n^ + BSyriy + CS'n, + H{t) + ....}] 
+ \ \^^^S + AS'^fix + BSffiy + CS^n^ + H{t) + ...)^ x S{r,t) , 

= S X i?e// (12) 

S{r,t)={S-{r,t),Sy{r,t),S\r,t)), S^ = l. (13) 

In fact, Eq.(12) was deduced from phenomenological grounds for bulk magnetic 
materials by Landau and Lifshitz originally in 1935. 

(c) Hamiltonian Structure of the LLG Equation in the Absence of Damping 

The dynamical equations for the lattice of spins (9) in one dimension (as well as 
in higher dimensions) posses a Hamiltonian structure in the absence of damping. 
Defining the spin Hamiltonian 

Hs = -J2 ^-^^+1 + ^(^i)' + ^i^i')' + + l^i^{t).Si) + (14) 



and the spin Poisson brackets between any two functions of spin A and B as 

{A,B}= ^ (15) 



a, 0,7=1 



Article submitted to Royal Society 



6 



Lakshmanan 



one can obtain the evolution equation (9) for A = from ^ = {Si, Hg}. Here Gai3j 
is the Levi-Civitta tensor. Similarly for the continuum case, one can define the spin 
Hamiltonian density 

Hs = lii^Sf + A{S^)^ + B{Syf + C{S'f + {H.S) + Hdemag + ] (16) 

along with the Poisson bracket relation 

{5«(r,t),5^(r',t')} =g„^^ ^^^(f - r',t - t'), (17) 

and deduce the spin field evolution equation (12) for A = 0. 
Defining the energy as 

E=^j (iV[(V5)2 + A(S'=^)2 + B{Syf + C(S'^)2 + H.S+ ....] (18) 

one can easily check that 




Then when A > 0, the system is dissipative, while for A = the system is conser- 
vative. 



3. Spin Torque Effect and the Generalized LLG Equation 

Consider the dynamics of spin in a nanoferromagnetic film under the action of a 
spin current (Stiles & Miltat, 2006; Bertotti et al. 2009). In recent times it has 
been realized that if the current is spin polarized, the transfer of a strong current 
across the film results in a transfer of spin angular momentum to the atoms of 
the film. This is called spin torque effect and forms one of the basic ideas of the 
emerging field of spintronics. The typical set up of the nanospin valve pillar consists 
of two ferromagnetic layers, one a long ferromagnetic pinned layer, and the second 
one is of a much smaller length, separated by a spacer conductor layer [6], all of 
which are nanosized. The pinned layer acts as a reservoir of spin polarized current 
which on passing through the conductor and on the ferromagnetic layer induces 
an effective torque on the spin magnetization in the ferromagnetic film, leading to 
rapid switching of the spin direction of the film. Interestingly, from a scmiclassical 
point of view, the spin transfer torque effect is described by a generalized version 
of the LLG equation (12), as shown by Berger (1996), and by Slonczewski (1996), 
in 1996. Its form reads 

— = Sx[H,ff + Sxj], S=iS-,Sy,S'), S^ = l (19) 
where the spin current term 

^- '^'^ np)-^-^. (20) 



f{P){3 + S.Sp) 4Pi 
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Here Sp is the pinned direction of the polarized spin current which is normally taken 
as perpendicular to the direction of flow of current, a is a constant factor related 
to the strength of the spin current and f{P) is the polarization factor deduced 
by Slonczewski (1996) from semiclassical arguments. From an experimental point 
of view valid for many ferromagneic materials it is argued that it is sufficient to 
approximate the spin current term simply as 

J = aSp (21) 

so that LLG equation for the spin torque effect can be effectively written down as 

(99 -. -. -, 

— =Sx[Heff + aSxSp], (22) 

where 

_ _ _ _ _ _ _ f)Q 

Heff = (V'5 + AS'^i + BS^j + CS'k + Hdemag + H{t) + ) - \S X —. (23) 

Note that in the present case, using the energy expression (18), one can prove that 

= j[-\\St\^ + a{StxS).Sp]d^r. (24) 



dE 
'dt 



This implies that energy is not necessarily decreasing along trajectories. Conse- 
quently, many interesting dynamical features of spin can be expected to arise in 
the presence of the spin current term. 

In order to realize these effects more transparently, let us rewrite the generalized 
LLG equation(22) in terms of the complex stereographic variable to (r, t) (Lakshmanan 
& Nakumara, 1984) as 

s- + isy 

- = TTT^' (25a) 
ax_ w + w* (^ ~ , 5^ _ (1 -^^*) ^25b) 



(l + ww*)' i{l+u)uj*)' (l+ww*) 

so that Eq.(20) can be rewritten (for simplicity Hdemag — ) 

2w*(Va;)2 A {1 - ij^){oj + CO*) 



i{l - iX)uJt + V a; - 



(l+ww*) 2 (1 + ww*) 



B {l + u:^){u:-u*) _ 1 - a;^;* 
2 (l + ww*) ^l+ujuj*' 

+ i(if- _ _ ^ 1 ^ ^■y^^■^ ^ ^2^ _ ^ ^ (26) 

where j = aSp, ujt = (ff )• 

It is clear from Eq.(26) that the effect of the spin current term j is simply to 
change the magnetic fleld H = {H'^,Hy,H'') as [H^ - if,Hy +ijy,H''+if). 
Consequently the effect of the spin current is effectively equivalent to a magnetic 
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field, though complex. The consequence is that the spin current can do the function 
of the magnetic field perhaps in a more efficient way because of the imaginary term. 

To see this in a simple situation, let us consider the case of a homogeneous 
ferromagnetic film so that there is no spatial variation and the anisotropy and 
demagnetizing fields are absent, that is we have (Murugesh & Lakshmanan, 2009a) 



(1 — iX)uJt = —{a — iH^)LO. 



(27) 



Then on integration one gets 
w(t) = ijj{Q)exp[ 



= u){0)exp 



-{a - iir-)t 

(1-*A) J 



exp 



.{a\-H^)t 



(1 + IAP) 



(28) 



(1 + IAP) 

Obviously the first exponent describes a relaxation or switching of the spin, while 
the second term describes a precession. From the first exponent in (28), it is clear 
that the time scale of switching is given by . Here A is small which implies 

that the spin torque term is more effective in switching the magnetization vector. 
Furthermore letting term to become zero, we note that in the presence of damp- 
ing term the spin transfer produces the dual effect of precession and dissipation. 

In Fig. 2, we point out clearly how the effect of spin current increases the rate of 
switching of the spin even in the presence of anisotropy. Further, one can show that 
interesting bifurcation sccnerio, including period doubling bifurcations to chaotic 
behaviour, occurs on using a periodically varying applied magnetic field in the 
presence of a constant magnetic field and constant spin current (Murugesh & Lak- 
shmanan, 2009a) . Though a periodically varying spin current can also lead to such 
a bifurcations-chaos scenerio (Yang et al. 2007), we believe the technique of apply- 
ing a periodic magnetic field in the presence of constant spin current is much more 
feasible experimentally. To realize this one can take (Murugesh & Lakshmanan, 
2009b) 



-ffe// = K(5.e||)e|| + Hdernag + H{t), 



(29) 



where k is the anisotropy parameter and e|| is the unit vector along the anisotropy 
axis and 



emag 



Choosing 



H(t) = {0,0, H^), 6*11 = {sin0iicos(j)ii, sinO\\sin(j)\\, cos0\i). 



(30) 



(31) 



the LLG equation in stereographic variable can be written down(in the absence of 
exchange term) as (Murugesh & Lakshmanan, 2009a,b) 



(1 — iX)uJt = — j{a — iH^)u) + iS'||K7 
z47r7 



(1 + M2) 



^1, 



Nail - \iof)uj - -^{1 -co^- Iw^o; 



-^(l+u;^-M> 



(32) 
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h„! (Oe) 

Fig.2:Switching of spin due to spin current in the presence of anisotropy 
(Murugesh & Lakshmanan, 2009a). 
Fig. 3: Bifurcation diagram corresponding to Eq.(32). Here ha^ = (Murugesh & 

Lakshmanan, 2009b). 

where 5*11 = S.e\\. Solving the above equation numericahy one can shown that 
Eq.(27) exhibits a typical period doubling bifurcation route to chaos as shown in 
Fig.3. 

The existence of periodic and chaotic spin oscillations in a homogeneous nano- 
spin transfer oscillator (STO) leads to exciting other possibilities. For example, 
an array/network of STO's can lead to the possibility of synchronized microwave 
power or synchronized chaotic oscillations (Grollier et al. 2006). Such studies are 
in progress. Other possibilities include the study of inhomogeneous films (including 
spatial variations) and discrete lattices, including higher dimensions (Bazaliy et al. 
2004). 

4. Anisotropic Heisenberg Spin Lattice 

Next we consider the dynamics of a discrete anisotropic Heisenberg spin systems 
without damping. Consider the Hamiltonian 

H = — ^(AS'j^S'i^^j + BS'j'^S^^_^_l + CS^S^^i) — Z? ^(S*^)^ — H. ^ Sn (33) 

n 

so that the equation of motion becomes (using the Poisson bracket relations (15)) 

— Sn'x[A{S''^j^i + S^_i)i + B{S'^^j^i + S^_i)j 

+ C(fi^+i + S^^_^)k + 2DS^k] + SnXH, (34) 
n = 1,2,3,...N. 

Recently (Lakshmanan & Saxena, 2008; Roberts & Thompson, 1988) it has been 
found that the coupled system (34) admits several classes of exact solutions, though 
the system may not be completely integrable for any choice of the parameters, 
including the pure isotropic one {A = B = C — 1,D — 0,H — 0). 
The class of exact solutions to (34) are as follows. 
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(a) Spatially homogeneous time dependent solution: 



SI = -y^l_72fc'2 sn{ujt + 5,k), (35a) 

SI = v'l -7^ cn{ujt + 5,k), (35b) 
5^ = 7 dn{ujt + S,k), (35c) 



where the frequency 



u> = 2'y^{B-C){A-C), k' = L_L^^_^ iB>A>C). (35d) 

Here 7 and (5 are arbitrary parameter and k is the modulus of the Jacobian eUiptic 
functions. 

(6) Spatially oscillatory time periodic solutions 



SI = (-1)"+Vl -7^fc'^ sn{ujt + 5,k), (36a) 
5^ = (-l)"v/l-7' cn{ujt + 5, k), (36b) 
S^ = -f dn{ojt + 6,k), (36c) 



where 



a; = 27V(^ + C)(B + C), fc^ ^ (36d) 

This solution corresponds to a nonlinear magnon. 

(c) Linear magnon solutions 
In the uniaxial anisotropic case A = B < C, the linear magnon solution is 

S^ = Vl -7^ sin{pn -iJt + 6), (37a) 
5^ = Vl -7^ cos(pn - + 5) , (37b) 
5^=7 (37c) 

with the dispersion relation 

(^ = 2^(0- Acqs p). (37d) 

(rf) Nonplanar static structures for XYZ and XYY models 
In this case we have the periodic structures 



S^ = Vl-72fc'2 sn{pn + 5,k), (38a) 
Sy = -72 cn(pn + 5, fc), (38b) 
5^ = 7 dn(pn + 5,A;), (38c) 
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where 

^' = IJ3|J' dn{p,k) = j. (38d) 

In the limiting case k = 1, one can obtain the locaHzed single soliton (solitary wave) 
solution 



S^=tanh{pn + S), SI = ^1-7^ sech{pn + 5), 8^= sech{pn + 5). (39) 

(e) Planar (XY) Case: 



Case{i) : = sn{pn + S,k), = cn{pn + S,k), = 0, (40) 
where dn{p, fc) = ^. In the limiting case, A; = 1, we have the solitary wave solution 
S^=tanh{pn + 6,k), = sech{pn + 5), S^ = 0. (41) 

Case{ii): = k sn{pn + S,k), = dn{pn + 6,k), 3^=0, (42) 
where cn{p, k) = ^. In the limit fc = 1, we have 

S^ = tanh{pn + S), 3^=0, = sech{pn + 6). (43) 

(/) Nonplanar XYY Structures: 

We have 

S!^ = cn{pn + 6,k), = ^ sn{pn + S,k), = - 7^ sn{pn + d,k), (44) 
where dn{p, fc) = ^. In the limiting case, we have the domain wall structure 

S!^ = sech{pn + 5), = ^ tanh{pn + 5), S"^ = V^l - 7^ tanh{pn + S). (45) 

In all the above cases one can evaluate the energies associated with the different 
structures and their linear stability properties. For details, one may refer to (Lak- 
shmanan & Saxena, 2008). 

{g) Solutions in the presence of onsite anisotropy and constant external magnetic 

field 

(i) Onsite anisotropy, D y^O, H = 0, A,B,C ^0: 

All the three types of solutions (35), (36) and (37) exist here also, except that 
the parameter C has to be replaced by (C — D) in each of these equations on their 
right hand sides. 
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(ii) Constant external field case, H = {H^, 0, 0), B = C ^ A, D = 0. 
An exact solution is 

= sn{pn + S,k), (46a) 
= sin{iut + 7) cn{pn + (5, fc), (46b) 
5^ — cos{u!t + 7) cn{pn + d, k), (46c) 

where dn{p, fc) = and a; = H^. 

One can study the linear stability of static solutions and investigate the existence 
of the so called Pcicrls-Nabarro barrier, that is whether the total lattice energy 
depends on the location of the soliton or not, for details see (Lakshmanan & Saxena, 
2008). 



{h) Integrability of the static case 

Granovskii and Zhedanov(1986) have shown that the static case of the pure 

anisotropic system 

Sn X [A(S^+i + S^_,) + B{Sl^, + Sl_,) + C{Sl^, + Sl_,)] = (47) 

is equivalent to a discretizcd version of the Schrodinger equation with two level 
Bargmann type potential or a discrete analog of Neumann system (Veselov, 1987) 
and is integrable. 



(z) Ishimori spin chain: 

There exists a mathematically interesting spin chain which is completely inte- 
grable and which was introduced by Ishimori (1982). Starting with a Hamiltonian 
H = —log{l + Sn-Sn+i), the spin equation becomes 



S„ — S„ X 



'n+l _|_ Sn- 



1 + 1 + 



(48) 



It admits a Lax pair and so is completely integrable. However, no other realistic spin 

system is known to be completely integrable. It is interesting to note that Eq.(48) 
also leads to an integrable reversible map. Assuming a simple time dependence, 

Sn{t) = {cos(j)ncosojt,cos^nsiniot,sin^n), (49) 

Quispel, Robert and Thompson (1988) have shown that Eq.(48) reduces to the 
integrable map 

Xn+l =[2x^ + iOxl + 2x„ - W - X„_i(-X^ - LOxl + LOXn + 1)] 

X [-x^-u)xl+ujXn + l-Xn-i{cox^-2xl-u;xl-2xn)]~^- (50) 

Finally, it is also of interest to note that one can prove the existence of localized 
excitations, using implicit function theorem, of tilted magnetization or discrete 
breathers (so called nonlinear localized modes) in a Heisenberg spin chain with easy- 
plane anisotropy (Zolotaryuk et al. 2001). It is obvious that there is much scope 
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for detailed study of the discrete spin system to understand magnetic properties, 
particularly by including Gilbert damping term and also the spin current, see for 
example a recent study on the existence of vortices and their switching of polarity 
on the application of spin current (Sheka et al. 2007). 

5. Continuum Spin Systems in (1+1) Dimensions 

The continuum case of the LLG equation is a fascinating nonlinear dynamical 

system. It has close connection with several intcgrablc soliton systems in the absence 
of damping in (1+1) dimensions and possesses interesting geometric connections. 
Then damping can be treated as a perturbation. In the (2+1) dimensional case novel 
structures like line soliton, instanton, dromion, spatiotemporal patterns, vortices, 
etc. can arise. They have both interesting mathematical and physical significance. 
We will briefly review some of these features and indicate a few of the challenging 
tasks needing attention. 

(a) Isotropic Heisenberg Spin System in (1+1) Dimensions 

Considering the (1+1) dimensional Heisenberg ferromagnetic spin system with 
nearest neighbour interaction, the spin evolution equation without damping can be 
written as (after suitable scaling) 

St = SxS,,, S = {S^,Sy,S'), S^ = l. (51) 

In Eq.(51) and in the following suffix stands for differentiation with respect to that 
variable. 

We now map the spin system (Lakshmanan et al. 1976) onto a space curve (in 
spin space) deflned by the Serret-Prenet equations, 

eix = D X Ci, D = Tei+ nes, Ci.ei = 1, i = 1, 2, 3, (52) 

where the triad of orthonormal unit vectors ei, 62, 63 are the tangent, normal and 
binormal vectors, respectively, and x is the arclength. Here n and r arc the curvature 
and torsion of the curve, respectively, so that = e\x-Six (energy density), k^t = 
ei.{eix X eixx) is the current density. 

Identifying the spin vector S{x,t) of Eq.(51) with the unit tangent vector ei, 
from (51) and (52), one can write down the evolution of the trihedral as 

eit = n X Cj, Ci = {ui, 0)2,1^3) = (^-^ -'r'^,-Kx,-HT^ ■ (53) 
Then the compatibility {ei)xt = {ei)tx, i = 1,2,3, leads to the evolution equation 

Kt = -2kxT - KTx, (54a) 
rt = {— -t'^)x + KKx, (54b) 

which can be rewritten equivalently (Lakshmanan, 1977) as the ubiquitous nonlin- 
ear Schrodinger (NLS) equation, 

iqt + qxx + 2\q\'^q = o, (55) 
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through the complex transformation 



1 

q = —K exp 



+00 



T 



dx 



(56) 



and thereby proving the complete integrability of Eq.(51). 

Zakharow and Takhtajan (1979) have also shown that this equivalence between 
the (1+1) dimensional isotropic spin chain and the NLS equation is a gauge equiv- 
alence. To realize this, one can write down the Lax representation of the isotropic 
system (51) as (Takhtajan, 1977) 

(l)ix = Ui{x,t,X)(j)i, (j)it = Vi{x,t,X)(j)i, (57) 

where the (2 x 2) matrices Ui = iXS, Vy = XSS^ + 2iX^S, S = (^^l ^^^^ , 
= 5^ ±15". Then considering the Lax representation of the NLS equation (55), 

(f>2x = U2(t)2, 4>2t = V2<^2, (58) 

where U2 = (^0 + AAi), V2 = {Bq + XBi + A^Bs), 

-^(-, '»•)• 

Bi = 2Ao, B2 = 2Ai, CTi's are the Pauli matrices, one can show that with 
the gauge transformation 

= 5"V2, S = (59) 

(57) follows from (58) and so the systems (51) and (55) are gauge equivalent. 
The one soliton solution of the component can be written down as 

S'^tx, t) = l- ^ „ sech'^^(x - 2r]t - x°), ^, rj, x° : constants (60) 

Similarly the other components and 5" can be written down and the N-soliton 
solution deduced. 

(6) Isotropic chain with Gilbert damping 
The LLG equation for the isotropic case is 

St = SxS,, + X[S,, - {S.S,,)§\. (61) 

The unit spin vector S{x, t) can be again mapped onto the unit tangent vector ei 
of the space curve and proceeding as before (Daniel & Lakshmanan, 1983) one can 
obtain the equivalent damped nonlinear Schrddinger equation, 

/ 

(99*' ~ ], (62) 

-00 

where again q is defined by Eq.(55) with curvature and torsion defined as before. 
Treating the damping terms proportional to A as a perturbation, one can analyse 
the effect of damping on the soliton structure. 
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(c) Inhomogeneous Heisenberg Spin System 

Considering the inhomogeneous spin system, corresponding to spatially depen- 
dent exchange interaction, 

St = (72 + fJ-2x)S X Sxx + A*2<5 'X Sx- (71 + y'ix)Sx, (63) 

where 71, 72, /Ui and 112 are constants, again using the space curve formalism, 
the present author and Robin BuUough (Lakshmanan & Robin Bullough, 1980) 
showed the geometrical/gauge equivalence of Eq.(63) with the linearly a;-dependent 
nonlocal NLS equation, 

iqt = ifxiq + i(7i + /Uia;)^^ 

\q\'^dx) = 0. (64) 

-00 

It was also shown in (Lakshmanan & Bullough, 1980) that the both the systems 
(63) and (64) are completely integrable and the eigenvalues of the associated linear 
problems are time dependent. 

(rf) n- dimensional Spherically Symmetric (radial) Spin System 

The spherically symmetric n -dimensional Heisenberg spin system (Daniel et al. 
1994) 



St{r,t) = Sx 



r 



(65) 



S2(r,i) = l, 5=(5^5^5"), =rf + r^ + ...+<, < r < 00, 



can be also mapped onto the space curve and can be shown to be equivalent to the 
generalized radial nonlinear Schrddinger equation, 

iqt + qrr + = - 2|?|' - 4(n - 1) Kdr) q. (66) 



It has been shown that only the cases n = 1 and n = 2 are completely integrable 
soliton systems (Mikhailov & Yaremchuk, 1982; Porsezian & Lakshmanan, 1991) 
with associated Lax pairs. 

(e) Anisotropic Heisenberg spin systems 

It is not only the isotropic spin system which is integrable, even certain anisotropic 
cases are integrable. Particularly the uniaxial anisotropic chain 

St = Sx [S^^ + 2 A S'n, + H], = (0, 0, 1) (67) 

is gauge equivalent to the NLS equation (Nakamura & Sasada, 1982) in the case 
of longitudinal field H = {0,0, H^) and is completely integrable. Similarly the 
bianisotropic system 

St = SxJSxx, J = diag{Ji,J2,Jz), J\i^Jii^Ji, (68) 
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possesses a Lax pair and is integrable (Sklyanin, 1979). On the other hand the 
anisotropic spin chain in the case of transverse magnetic field, H = {H^,0,0), 
is nonintegrable and can exhibit spatiotemporal chaotic structures (Daniel et al. 
1992). 

Apart from the above spin systems in (1+1) dimensions, there exists a few other 
interesting cases which are also completely integrable. For example, the isotropic 
biquadratic Heisenberg spin system, 



Sx 



^xx ~i~ T 1 ^XXXX "^{^ -^XX^^XX "^{^ '^XXX^^X 



(69) 



is an integrable soliton system (Porsezian et al. 1992) and is equivalent to a fourth 
order generalized nonlinear Schrddinger equation, 

iqt + Qxx + 2|5r + 7 [Qxxxx + 8\q\^qxx + 2q^q*^ + 4g|feP + &q*ql + 6\q\^q] = 0. 

(70) 

Similarly the SO{3) invariant deformed Heisenberg spin equation, 

St = SxS,, + aS,{S,f, (71) 

is geometrically and gauge equivalent to a derivative NLS equation (Porsezian et al. 
1987), 

iqt + qxx + ^\q\'^q - ia{\q\'^q)x = 0. (72) 

There also exist several studies which maps the LLG equation in different limits 
to sine-Gordon equation (planar system), Korteweg-de Vries (KdV), mKdV and 
other equations, depending upon the nature of the interactions. For details see for 
example (Mikeska & Steiner, 1991; Daniel & Kavitha, 2002). 

6. Continuum Spin Systems in Higher Dimensions 

The LLG equation in higher spatial dimensions, though physically most important, 
is mathematically highly challenging. Unlike the (1+1) dimensional case, even in 
the absence of damping, very few exact results are available in (2+1) or (3+1) 
dimensions. We briefly point out the progress and challenges. 

(a) Nonintegrability of the Isotropic Heisenberg Spin Systems 
The (2+1) dimensional isotropic spin system 

St = Sx (^S^^ + Syy), S={S^,Sy,S'), 52 = 1, (73) 

under stereographic projection, see Eq.(25), becomes (Lakshmanan & Daniel, 1981) 

(1 + ujoj*)[ioJt + Wxx + i^yy] - 2w*(w^ + w^) = 0. (74) 
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It has been shown (Senthilkumar et al. 2006) to be of non-Painleve nature. The 
solutions admit logarithmic type singular manifolds and so the system (73) is non- 

integrable. It can admit special types of spin wave solutions, plane solitons, ax- 
isymmetric solutions, etc. (Lakshmanan & Daniel, 1981). Interestingly, the static 
case 

2u>* 

<^xx + Wto = (i^^^^*)2 (^x + ^l) (75) 
admits instanton solutions of the form 

with a finite energy (Belavin & Polyakov, 1975; Daniel & Lakshmanan, 1983). Fi- 
nally, very little information is available to date on the (3+1) dimensional isotropic 
spin systems (Guo & Ding, 2008) 

(6) Integrable (2+1) Dimensional Spin Models 

While the LLG equation even in the isotropic case is nonintcgrablc in higher 
dimensions, there exists a couple of integrable spin models of generalized LLG 
equation without damping in (2+1) dimensions. These include the Ishimori equa- 
tion (Ishimori, 1984) and Myrzakulov equation (Lakshmanan et al. 1998), where 
interaction with an additional scalar field is included. 

(i) Ishimori equation: 



St — S y. {Sxx "I" ^yy) ~^ UxSy + UySx, (77a) 
Uxx - o-'^Uyy = -2a^S.{Sx X Sy), = ±1. (77b) 

Eq.(77) admits a Lax pair and is solvable by inverse scattering transform (d-bar) 
method (Konopclchcnko & Matkarimov. 1989). It is geometrically and gauge equiv- 
alent to Davey-Stewartson equation and admits exponentially localized dromion so- 
lutions, besides the line solitons and algebraically decaying lump soliton solutions. 
It is interesting to note that here one can map the spin onto a moving surface 
instead of a moving curve (Lakshmanan et al. 1998). 

(ii) Myszakulov equation I (Lakshmanan et al. 1998) 
The modified spin equation 

St = (SxSy + uS) ^ , Ux = -S.{SxxSy) (78) 

can be shown to be geometrically and gauge equivalent to the Calogero-Zakharov- 
Strachan equation 

iqt = Qxy + Vq, Vx = 2{\q\% (79) 
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and is integrable. It again admits line solitons, dromions and lumps. 

However from a physical point of view it will be extremely valuable if exact 
analytic structures of the LLG equation in higher spatial dimensions are obtained 
and the so called wave collapse problem (Sulem & Sulem, 1999) is investigated in 
fuller detail. 



(c) Spin Wave Instabilities and Spatio-temporal Patterns 

As pointed out in the beginning, the nonlinear dynamics underlying the evolu- 
tion of nanoscale ferromagnets is essentially described by the LLG equation. Con- 
sidering a 2D nanoscale ferromagnetic film with uniaxial anisotropy in the presence 
of perpendicular pumping, the LLG equation can be written in the form (Kosaka 
et al. 2005) 

St = Sx F,ff -XSx—, (80a) 

where 

Feff = J^^S + Ba + KS\\e\\+Hm, (80b) 

Ba = ha±icosujt 1 + sinLot j) + /lay 6*11 (80c) 

Here i, j are unit orthonormal vectors in the plane transverse to the anisotropy 
axis in the direction ey = (0, 0, 1), k is the anisotropy parameter, J is the exchange 
parameter, Hm is the demagnetizing field. Again rewriting in stereographic coordi- 
nates, Eq.(80) can be rewritten (Kosaka et al. 2005) as 

i(l - iX)u}t = J V^w - - ha\\ -v + iav + K- ' ' ' 



(l+ww*); V (l + l'^l') 

+\k^{1 - - K\\u^ + \ {Hme-'"' - v^H*^e'-') co. (81) 

Then four explicit physically important fixed points (equatorial and related 
ones) of the spin vector in the plane transverse to the anisotropy axis when the 
pumping frequency v coincides with the amplitude of the static parallel field can be 
identified. Analyzing the linear stability of these novel fixed points under homoge- 
neous spin wave perturbations, one can obtain a generalized Suhl's instability cri- 
terion, giving the condition for exponential growth of P-modes (fixed points) under 
spin wave perturbations. One can also study the onset of different spatiotemporal 
magnetic patterns therefrom. These results differ qualitatively from conventional 
ferromagnetic resonance near thermal equilibrium and are amenable to experimen- 
tal tests. It is clear that much work remains to be done along these lines. 



7. Conclusions 

In this article while trying to provide a bird's eye view on the rather large world of 
LLG equation, the main aim was to provide a glimpse of why it is fascinating both 
from physical as well as mathematical points of view. It should be clear that the 
challenges are many and it will be highly rewarding to pursue them. What is known 
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at present about different aspects of LLG equation is barely minimal, whether it 
is the single spin case, or the discrete lattice case or the continuum limit cases 
even in one spatial dimension, while little progress has been made in higher spatial 
dimensions. But even in those special cases where exact or approximate solutions are 
known, the LLG equation exhibits a very rich variety of nonlinear structures: fixed 
points, spin waves, solitary waves, solitons, dromions, vortices, bifurcations, chaos, 
instabilities and spatiotemporal patterns, etc. Applications are many, starting from 
the standard magnetic properties including hysterisis, resonances, structure factors 
to applications in nanofcrromagncts, magnetic films and spintronics. But at every 
stage the understanding is quite imcomplete, whether it is bifurcations and routes 
to chaos in single spin LLG equation with different interactions or coupled spin 
dynamics or spin lattices of different types or continuum system in different spatial 
dimensions, excluding or including damping. Combined analytical and numerical 
works can bring out a variety of new information with great potential applications. 
Robin BuUough should be extremely pleased to see such advances in the topic. 
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